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We propose a simple ansatz that allows to generate new exactly solvable multi-state Landau- 
Zener models. It is based on a system of two decoupled two-level atoms whose levels vary with time 
and cross at some moments. Then we consider multiparticle systems with Heisenberg equations for 
annihilation operators having similar structure with Shrodinger equation for amplitudes in multistate 
Landau-Zener models and show that the corresponding Shrodinger equation in multiparticle sector 
belongs to the multistate Landau-Zener class. This observation allows to generate new exactly 
solvable models from already known ones. We discuss possible applications of the new solutions in 
the problem of the driven charge transport in quantum dots. 



Landau-Zener (LZ) theory [1,2] is one of the most important and influential results in non-stationary quantum 
mechanics. Last decade a generalization of LZ theory to more than two states attracted particular attention due to 
numerous applications in atomic and molecular physics [3], [4], nanomagnets [5], Bose-Einstein condensate [6] and 
systems with avoided band crossings [7], [8]. The multi-state LZ problem (see for example [9] ) is concerned with 
finding of the transition amplitudes for a system with the Hamiltonian, whose matrix form reads: 



where _B is a diagonal matrix and the matrices A and B are independent of time. In its general form this problem is 
still unsolved, but a number of exact results for special choices of the matrices B and A were found [9-19]. 

In almost all available exact solutions the transition probabilities are expressed in terms of the genuine two-level 
LZ formula successively applied at each diabatic level intersection. In other physical problems such a procedure is 
often applied as an approximation. These problems include atomic and molecular collisions [20] and the transitions 
at crossing of two Rydberg multiplets of energy levels [3] . 

In this work we find very simple ansatz that generates new solvable models and may explain the properties of 
already known solutions. The main idea employs the single-particle Hamiltonian which acts independently in several 
two-dimensional subspaces of the Hilbert space. It is worth mentioning that while results in one-particle sector are 
trivial, the same Hamiltonian generates non-trivial solutions in the many particle spaces. Such a construction is 
akin to the group-theoretical method of finding higher irreducible representations as a symetrized direct product of 
the fundamental representation. Using this method we can study the problem of driven charge transport through a 
quantum dot and find new solutions in multistate LZ theory. Particularly, we derive the transition probabilities for a 
four state LZ problem which is very similar to the four state bow-tie model and for a problem of intersection of two 
bands of parallel levels. 

This article is organized as follows: in section II we show how already known solutions of LZ models can generate new 
exactly solvable models with the Hamiltonian (1). We demonstrate how the exact solution for two independent two 
level systems can generate a new solution of a four-state LZ model. In III we generalize Demkov-Osherov solution to 
the case of many particles and use the result for derivation of master equations that describe a driven charge transport 
in quantum dots. In IV we provide an example of a solvable model that can be generated from the Demkov-Osherov 
solution. 



I. INTRODUCTION. 



H = Bt + A, 



(1) 



II. BOSONIC MULTI-STATE LZ MODELS. 



Lets consider a Hamiltonian that describes the interaction of four bosonic fields a, b, c, d: 



H = {f3it + Ei)a+a + (pst + Ej,)d+d+ (/Ja* + £2)0+ c + {(3^1 + Ei)b+b+ 
g(a+b + b+d) + j(c+d + d+c) 



(2) 
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This Hamiltonian depends explicitly on time and conserves the total number of particles in the system. Therefore it 
can be considered independently in subspaces with fixed total number of particles. Let |0 > be the vacuum state. The 
Hamiltonian (2) describes the evolution of two disjointed systems. However, being projected onto the 2-particle sector, 
its matrix form looks less trivial. The complete two-particle sector is the 10-dimcnsional Hilbert space spanned onto 
direct products of any two single-particle states. The four-dimensional subspace R4 of the 2-particle sector spanned 
onto vectors: 



|2 >= a+d+\0 > 
\3>=d+b+\0> 
14 >= c+S+IO > 



(3) 



is invariant with respect to the action of the Hamiltonian (2). Hence, if the initial state belongs to this subspace, the 
state vector at any time remains in R4: 



\m) - c,{t)\l) + C2{t)\2) + C3WI3) + Ciit)\A) 
In the basis (3) the Hamiltonian (2) has the following 4x4 matrix form; 
/ iPi + P2)t + {El + E2) 7 



H = 



7 


\9 



9 

iPi + p3)t + {El + E3) g 
g {Pa + 0i)t + {E3 + E4) 7 

7 {P2+ I3i)t + {Ei + E2) J 



(4) 



(5) 



The problem described by the Hamiltonian (5) belongs to the multistate Landau-Zener class (1). 

We should point out that it cannot be mapped on the already known exactly solvable multistate LZ models. 




FIG. 1. Time dependence of the adiabatic energies (solid lines) and diagonal elements (dashed lines) of the Hamiltonian (5). 
The choice of parameters is /3i = 5, (32 = —3, jSs = 0, /34 = —1.5, Ei = 3, E2 = 0.5, E3 = —2, E4 = —1.5, g = 1, 7 = 1.5 

In Fig.l we show the adiabatic energies of the Hamiltonian (5) as functions of time for typical choice of parameters. 
It is clearly seen that there is a point of an exact adiabatic level crossing (diabolic point) on the figure. 

In the Heisenberg representation the evolution equations decouple into two pairs of equations for bosonic operators: 

ih = {Pit + Ei)a + gb 

A . (6) 

ib = {Pit + E4)b + ga 

and 

i£ ^ {P2t + E2)c + jd 
id = {Pst + E3)d + 
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Let do, &07 Co, do denote the operators a, b, c, d at the initial moment of evolution. Then the solutions of equations (6) 
and (7) are: 

d{t) = Sii{t)ao + Si2{t)bQ , , 

b{t) = S2l{t)ao + S22{t)bo 



cit) = S[^(t)co + S[2{t)do 

d{t) = S'2,{t)C0 + S'22it)do ^' 

Here Sij and are the matrix elements of the evolution operators for (6) and (7), respectively. Due to the linearity 
they arc the same for the opc;rator and numerical functions obeying these differential equations. Hc;nc;c;, we can extract 
them directly from the solution of the two-state LZ problem. For the evolution from t = —00 to t = +00 their squares 
of modulus are: 



Pi = 






l'S'22 




= e 


-2V/l/3i 


-M 


qi = 


\Sl2 


2 _ 


1^21 


2 


= 1 


-Pl 




P2 = 






\^22 




= e 


-27r7V|/32 


-/33I 


92 = 


1^12 








= 1 


-P2 





(10) 



Returning to the four-state LZ problem in the two-particle sector considered earlier, we first note that each state I7) of 

this subspace is the direct product of states from two independent subspaces of the one-particle sector \j) = \aj) 

aj = 1,2; =3,4 (note that here 1,2,3,4 enumerate single-particle state, for example |1) = a+|0)). The evolution 

matrix is also the direct product of evolution matrices in the independent subspaces of the one-particle sectors: 

U{t) = Ua{t) ^Un{t). Therefore transition matrix elements and probabilities Pij in the considered subspace are 

factorized: 



Pij — PdiajPfiifij (11) 

In terms of the LZ probabilities for two-level problems introduced earlier the transition probability matrix P, whose 
elements are defined by equation (11), reads: 

/ P1P2 pi92 gig2 pi92 \ 
p ^ Pi92 P1P2 q\P2 qiq2 ^2) 
9192 P2qi P1P2 piq2 

\ qiP2 9i92 ^192 P1P2 J 

This result does not depend on the parameters Ei. It is interesting that scattering matrices Sij {t) and 5^ {t) are known 
for any t [2] which make it possible to find the evolution operator at any time in the Schrodinger representation. 



III. DRIVEN CHARGE TRANSPORT THROUGH QUANTUM DOTS. 



In similar fashion to the previous section the fermionic systems can lead to Heisenberg equations for annihilation 
operators that have the same structure as Shrodinger equation for amplitudes for some exactly solvable multistate 
Landau-Zener model. As we will show, in the Shrodinger rc;presentation such a Fermi system with fixed number of 
particle is equivalent to a new solvable multistate Landau-Zener model. The models that we will examine correspond 
to the driven charge transport in nanostructures. 

Consider a quantum dot coupled to an external reservoir like the system shown in Fig2. Lets consider that initially 
some of the reservoir energy levels are filled with electrons, the others are empty. Lets assume the dot has only one 
electron bound state whose energy in real semiconductors can be regulated by the gate voltage; therefore the variation 
of the gate voltage with time generates time dependence of the dot's electronic level. 
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FIG. 2. A single energy level in a potential well coupled to two leads at zero temperature. Electron states in leads are filled 
up to Fermi energies, that can be difi^erent in right and left leads. 



The Hamiltonian of the electrons in the dot and reservoir reads: 



N 



(13) 



n=l 



Here cq is the fermionic operator that annihilates the electron on the dot level and c„ is the annihilation operators 
for the level En of the reservoir; E{t) is the time-dependent energy of the dot state. In our treatment the last term 
in (13) describes the tunneling between the leads and the single level in the quantum dot. We ignore all interactions 
among electrons except the one due to Pauli principle. 

Similar time-dependent single-particle problems for quantum dots have been already considered in [21]. Though 
our system is simplified but rather it is interesting because it has an exact solution. 

In the context of LZ theory, we approximate the dot energy by a linear function of time: E{t) = (3t. The Heisenberg 
operator equations corresponding to Hamiltonian (13) are: 



iCo = 0tCo + QnCn 
n 

iCn = EnCn+ QnCQ 

Due to the linear structure of these equations the solution can be formally written in the matrix form: 

d{t) = S{t)d{to) 



(14) 



(15) 



where c= (co,ci,...,c„) 

As in the previous section, the evolution matrix S{t) is completely determined by the coefficients of the differential 
equations (14) and is the same for operator and c-function solutions. Hence, it is enough to solve (14) with all operators 
replaced by c-functions. Such a system of equations coincides with that of the Demkov-Osherov model [18]. The latter 
provides transition amplitudes for a single energy level crossing an energy band consisting of time-independent levels. 

In Demkov-Osherov model the Shrodingcr equation for the amplitudes of different quantum states can be written 
as follows: 



/ «o(i) \ 
ai{t) 



( Pt 71 
7i* Oil 



In \ 



an / 



( ao{t) \ 



(16) 



where ak are ordered as follows ai < a2 < ■■■ < ctJVa (assuming that none two of the are equal; we also assume for 
definiteness that (3 > 0). 



The absolute values of the /S-matrix components {Ski = \c-k\'^)\/ O'l'^' [18] are: 
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S'oo = e-'^(^i+-+^") 

Soi = (1 - e-2'^^')i/2e-'r(^i+i+-+^n) 

Sko = e— - e-2--'=) 1/2, (fc = 1, . . . , n) . . 

S'fe/ = 0, (1 < A; < Z) ^ 

Sii = e-^^' 

Ski = (1 - e-2'^^')^^^e-'^(^'+i+-+^'=-i)(l - e-2'^^'')V2, (fc > /) 
(where the index I = 1. . .n and Zk = |7feP//3) 

The probabiUties to find an electron on a particular n-th level are. 

Pn =< C+{t +00)c„(i ^ +00) >= En2 Sln.Snn, < C+ (i ^ -Oo)c„,(i ^ -Oo) > = 

where 5^ = Sij(t — > +cxd) and the STinimation is taken over the initially filled states only. The scattering matrix 
elements S^^uf are given in (17). If the band of electron states in the external system is continuous then it is reasonable 
to use the approximation, in which g{E) = ^ while the value T{E) = 2wp{E)\g{E)\'^ is kept finite. Here p{E) is 

2 -En 

the density of states in the band and the elements of scattering matrix become \Soi\'^ = exp J ~^g^^ dE 

Now lets consider a dot that is connected to two leads. The left lead is characterized by the coupling function 

gL{E) and the densities p{^{E), p%{E) where f and e refer to the filled and empty states in the left lead (pl{E) — 
p{^{E) + pKE)), analogously wc can define the quantities guiE), p^j^{E), p'ji{E) for the right lead. Moreover it 
is more convenient to introduce the following notations: r{(i5) = 2tt p^j^{E)\gL{E)\^ , r|_(-E') = "^t^ P'L{E)\g l{E)\^ , 
4(£;) = 2np^^{E)\gn{E)\\V%{E) = 2np%{E)\gn{E)\\ = ^ + r{, and = V% + VI. 

If the dot state was initially empty and if this state crosses the region from energy Ei to energy £'2, then in 
the continuous approximation (18) leads to the following probability for the dot level to be finally filled after all 
Landau-Zener transitions: 

E2 

1 

'1. 



E2 , 

Pf{E2)=Po= / ^5-^ exp 



/3 



{rf{E)+r''{E))dE 



dE' (19) 



If the dot level was initially filled, it is necessary to add l^ooP = e ^1 to (19). One can check that 

the result (19) is the solution of the following system of differential equations: 

= -^%E)pf{E) + Tf{E)p,{E) 
P^ = -'^HE)p,{E)+T-{E)pf{E) ^ ' 

here Pe{E) = 1 — Pf{E) is the probability that the dot level will be empty when it has energy E. The equation for 
the charge that is transferred to the right lead can be derived in a similar way 

= {e/P){TUE)pAE) - TiiE)pAE)) (21) 

Note that equations (20), (21) were derived from the exact solution of the problem with microscopic Hamiltonian (13) 
rather than from random phase approximation or other type of phenomenology. 

Let us calculate the total charge transferred through the dot from the left lead to the right lead at zero temperature 
and a fixed bias that leads to a difference of Fermi energies in the left and in the right leads. Lets assume that the dot 
level was initially much lower than both Fermi levels and it was filled. Then the energy of this state grows linearly 
with time crossing both Fermi levels during the evolution. Since transitions will proceed presumably when the dot 
level is between Fermi energies of the leads, we can apply the following approximations: T'^j^{E) = Tl^I — 6{E — Cp)), 
THE) = Tl0{E - e^), Tf^{E) = rR{l - e{E - )) and T%{E) = TRe{E - e^) with Tr and Tl are constant. To 
find the total charge that is transferred to the right lead we formally put the final dot state energy equal to infinity 
in the solution of the equations (20) and (21). In the result the total charge transfered to the right lead is 

.(Ffl + Fi)^ f3 >^{Tr + Tl)\ ^ ' 

Clearly at Cp = Cp we find Q = erR/{rR + Fi,), which can be interpreted as the electron charge e multiplied by the 
probability for the electron that is initially placed into the dot to transfer to the right lead. 
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IV. SOLVABLE MODEL OF BANDS CROSSING. 



In this section we will construct a new solvable LZ model employing the fermionic Hamiltonian. The Hamiltonian 

(13) projected onto the /c-particle sector generates the evolution in the Hilbert space of dimensionality {N+l)l/{k\{N+ 
1 — fc)!). If we assume that the single-particle Hamiltonian laying in the background is the same as that of the Demkov- 
Osherov model, then all such models are reducible to this single-particle one. 

Generalized Landau- Zcncr models that deal with intersections of bands of parallel levels are important in many 
applications such as in driven tunneling of nanomagnets coupled to nuclear spins [5] and in driven charge transport 
in quantum dots [19]. 

Up to now only two exact solutions of this type were known: Demkov-Osherov solution and the case of the infinite 
number of states in bands that equally interact with states of another band [22], [10]. For an important case of a 
finite number of states in bands that is not equal to unity exact solutions for all transition probabilities have not been 
found yet though the absence of counterintuitive transitions was analytically proved [19]. Nearly-exact solution valid 
in the quasidegeneracy approximation was found and investigated in [15]. We will show that our method can be used 
to generate exactly solvable models with interband transitions. 

Lets consider a system of two Fermi particles with the Hamiltonian 



H = Eib+b + E2C+C + td+d + gi{d+d + d+a) + g2{b+d + d+b) + gsic+d + d+c) 



(23) 



Let > i^i > 0. As we demonstrated previously, the solution of the operator evolution equation can be written in 
the form: 



/ m \ 

ait) 

m 
V m J 



Sit,to) 



I dito) \ 

a(io) 



(24) 



where S{t,to) is the matrix of evolution for a 4-state Demkov-Osherov model. Lets restrict the Hilbert space to 
the subspace of only two particles. It includes six states: |1 >= rf+a+|0 >, |2 >= (i+6+|0 >, |3 >= rf+c+|0 >, 
|4 >= a+fo+lO >, |5 >™ a+c+lO > and |6 >= >. Similarly to the bosonic case, this subspace is invariant 

during the evolution process. The Hamiltonian restricted to this subspace has the following matrix form: 



H 











-92 


-93 





\ 





t + Ei 





91 





-53 










t + E2 





91 


92 




-92 


91 





El 










-93 





91 





E2 







Vo 


-93 


92 








El + E2 


J 



(25) 




FIG. 3. Diagonal elements of the Hamiltonian (25) as functions of time. 



Let Pij = 1, ...,6) be the probability to transit from the state j to the state i after the band crossing. The 
transition probabilities can be expressed in terms of the fermi-operators in the Heisenberg representation at f — > 00. 
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Pin =< n\a'^ ad^ d\n > 
P2„ =< n\b+bd+d\n > 
Psn =< n|c+crf+d|n > 

Pau —< n\d'^ab'^b\n > 
Pbn =< n\a^ac^c\n > 
Pg„ =< n\b+bc+c\n > 



(26) 



Substituting (24) into (26) and employing the elements of the evolution matrix from (17) we get the following result: 



P = 



where 



/ P2P3 


qiq2P3 


qm 


P192P3 Pi 93 \ 





PiPa 


Pl9293 


9iP3 919293 P293 








P1P2 


qiP2 92 


92 


qiP2 





P1P2 


P2q3 


919293 


qips 


P19293 P1P3 




pm 


92P3P1 


91 93 9192P3 P2P3 J 




Pi = e 








qi = '^ 


-Pi, {i 


= 1,2,3) 



(27) 



(28) 



V. CONCLUSIONS. 



In conclusion, we presented the procedure that generates new exactly solvable multi-state LZ models. Some of 
them are useful for the description of driven charge transport in quantum dots and driven tunneling in nanomagnets. 
As an example, we derived two new solvable models and found the transition probability matrices for them. 

There have been three known classes of solvable multi-state LZ models that provide transition probabilities for a 
finite number of states: 

1. The Demkov-Oshcrov model. 

2. The SU(2) symmetry class that deals with an arbitrary spin in external magnetic field with the following 
Hamiltonian: 

H = tS,+ gS^ (29) 



3. The generalized bow-tie model that treats the case when two levels are parallel while the other levels intersect 

at one point between the parallel ones. 

This list can be extended with different generalizations of these models to the case of degenerate states. For example, 
it is possible to solve the LZ model for two degenerate levels by changing basis in such a way that all equations 
decouple into independent two state Landau-Zencr transitions. It is worth mentioning that sometimes a few elements 
of transition probability matrix can be found while the others remain unknown [9], [19]. 

All these models provide very simple results. For example, transition probabilities in the Demkov-Osherov model 
coincide with those taken from successive application of the two state Landau-Zencr formula. The same is true for 
the generalized bow-tie model. Finally in all models the transition probabilities are simple polynomials of Zh = 
exp(— 7r|g(fc|^). This fact gives a strong feeling that there should be a common symmetry in the background of all these 
models. Our results demonstrate the same properties and we know that the reason for this was the symmetry that 
makes the Hamiltonian equivalent in some sense to the one for a much simpler problem. 

We did not study deeply the relations between the known models and our new solutions but there are indications 
that such a relation can exist. For example the model (5) and the generalized bow-tie model arc very similar. At 
g = J and /3i = —P2, P3 = —Pi the Hamiltonian (5) belongs to the class of generalized bow-tie models. Also we note 
that models of SU{2) class [13] can be derived from the following bosonic Hamiltonian: 

H = ta+a-th+b + g{d+b + b+d) (30) 
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In the singlc-particlc sector the Hamihonian (30) leads to the simple two-state LZ model. In the iV-boson sector the 
Schrodinger equation for diabatic states coincides with ones for a spin S = N/2 in magnetic fields. This construction 
is an application of the Schwinger bosons [23] to the LZ problem. 

It is interesting to check what models can be reduced to decoupled two level systems. Probably this can be done 
using of the group representation theory. 
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